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I. INTRODUCTION 

In mathematical sciences some equations are highly privileged in the sense that their incarnations 
in analogous forms [1], at times of course with duifferent meanings of the underlying symbols, explain 
altogether different phenomena in Nature. One such equation, besides the equation of continuity, 
is the diffusion-reaction (D-R) equation which has offered explanation of many phenomena lying 
in the domains of not only physics and chemistry but also biology and (now) perhaps social and 
economic sciences. While the linear version of the D-R equation under its various names such as heat 
equation, Fokker-Planck equation, Schrodinger-like equation, etc. has been studied very extensively 
in different contexts, the study of its nonlinear versions and sometimes in higher dimensions have 
also evoked considerable interest in recent years [2-8]. 

Before proceeding further some important remarks about the mathematical contents of the D-R 
and Schrodinger equations are worth making. The linear version of the D-R equation, namely 

- DV 2 C(x, t) + v.VC(x, t) + U(x)C(x, t) = y-L, (1) 

can be compared with the time-dependent Schrodinger equation 

^V 2 ^(x, t) + V(xMx, t) = i^^, (2) 

where h=m= 1. In eq.(l) D is the diffusion coefficient and the velocity v in general is a function 
of x and t. In the prsent work, however, we shall consider v as a constant, independent of both 
space and time. Note that in mathematical literature eqs. (1) and (2) are classified differently 
in view of the different long-time behaviour of their solutions. As a matter of fact, in the limit 
t — > oo, while the complete solution of (1) vanishes, the solution of (2) remains periodic in time. 
The presence of the velocity-term in (1) further makes it distinct from (2). This term, in fact, makes 
the corresponding 'Hamiltonian' non-Hermitian and as a result the problem of reality of eigenvalues 
of such a Hamiltonian becomes of considerable interest in view of some recent studies[9-ll] in this 
direction. We shall address some of these issues in the present work. Lastly, while eq.(2) in quantum 
mechanics is set on the basis of some physical requirements with a rich physics content in it, eq.(l), 
on the other hand, is just a classical one like any other partial differential equation in mathematics, 
of course with different contextual meanings of various symbols. 
In the present work, we shall investigate the solutions of one-dimensional linear D-R equation, 

C t + vC x = DC xx -V(x)C, (3) 
and some of its nonlinear versions, namely 

C t + vC x = DC XX + [a + U(x)]C - b\C\C, (4) 

and 

C t + vC x = DC XX + [a + U{x)]C - b\C\ 2 C, (5) 

where a and b are real constants. We were motivated to study these equations from the recent works 
of Moiseyev and Gluck[5] and of Nelson and Shnerb[3]. These authors study the three-dimensional 
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version of these equations with reference to the derealization problem in population biology- a fea- 
ture of the non- Hermit ian character of the Hamiltonian associated with these equations. No doubt, 
the present study of one-dimensional systems will have limited scope as far as the applications to 
physical problems are concerned, nevertheless it will provide some clue towards a better understand- 
ing of these 'real-world' problems in mathematical terms. While the eq. (4) can be considered as a 
particular type of generalization of the Malthus-Verhulst growth model [6] in the study of biological 
systems, eq.(5) is the analogue of nonlinear Schrodinger equation. These cases, to the best of our 
knowledge, have not been studied earlier. In fact, eq.(4) is a slightly different version of the equation 
studied by Nelson and Shnerb[3] in the sense that the quadratic nonlinearity now appears as \C\C 
in stead of C 2 . 
With regard to eq.(3), we recast it in the form 



HC{x,t) = — — , (6) 

where H, a non-Hermitian 'Hamiltonian' operator, is given by [5] 

H = -D— + v— + V(x). 
ox z ox 

To be more specific, one writes the solution of (3) as C(x,t) = ip(x).r(t) and uses the method of 
separation of variables to obtain the eigenvalue equations, 



Hi/>(x) = \i>(x), (7) 

with 



dr ■> 

solution with time. Note that for a complex eigenvalue A there is a possibility of retaining the 

periodic behaviour of solutions with time- a feature built already in the Schrodinger eq.(2). Further, 

A." 

d,v ■ 



and (^) = —At. The solution of the latter equation implies an exponential decay of the complete 
with time. Note that for a complex 
behaviour of solutions with time- a fea 
we write the 'classical' analogue (using p = —i4-) of H as 



H(x,p) = Dp 2 + ivp + V(x). (9) 

Note that for a symmetric potential the non-Hermitian Hamiltonian (9) is symmetric under parity 
(x — * —x) operation (P) and it is not symmetric under time-reversal (t — > —t, i — > —i) operation (T) 
or under the combined operations of parity and time reversal (PT). Therefore, in the light of the 
conjecture/prescription suggested by Bender et al[9] and used by others[10] eq. (7) can not admit 
real eigenvalues, as H. is not a PT-symmetric Hamiltonian though it is non-Hermitian. For the case 
of a symmetric potential V(x), we shall investigate this problem of reality of eigenvalues here within 
the framework of a complex phase space approach proposed[12] and used recently[13] for a variety 
of complex potentials. In this approach, the real (x,p)phase plane is extended to a complex phase 
space characterized by writing x and p as 

x = xi + ip 2 ;p = pi + ix 2 , (10) 
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where (xi,p±) and (X2, P2) turn[15] out to be the canonical pairs in an equivalent real two-dimensional 
space. The arrangement of the paper is as follows: 

In Sect. 2 we investgate the solitary wave solutions of the nonlinear eqs.(4) and (5). We study the 
eigenvalue problem associated with the non-Hermitian operator (8) in Sect. 3. In Sect. 4, we highlight 
some other, so far unexplored, mathematical features of the Hamiltonian (8) (or (9)). In particular, 
for the linear case, in analogy with the Schrodinger equation[16] an Ermakov system[17] of equations 
is derived in this section. Also, corresponding to the one-dimensional complex Hamiltonian (9), an 
equivalent system of two, two-dimensional real Hamiltonians is obtained for an analytic potential 
function V(x). Finally, concluding remarks are made in Sect. 5. 

II. SOLITARY WAVE SOLUTIONS OF THE NONLINEAR D-R EQUATIONS 

In this section we obtain solitary wave solutions of the nonlinear eqs.(4) and (5). 

A. Solution of eq.(4) 

For the traveling wave solution of eq.(4) we consider the case when the random potential U(x) is 
constant, say Uq, and for the solution C(x,t) we make an ansatz 

C(x,t)=r(O.exp\i0(Z) + 6t], (11) 

where £ = x — wt, and 5 and w are arbitrary real constants to be determined later. Using (11) in 
(4) and then separating the real and imaginary parts of the resultant expression, one obtains 

(v - w)r' + r5 = Dr" - Dr9' 2 + (a + U )r - br 2 , (12) 

(y -w)r9' = D(2r'6' + r9"). (13) 

Here, the primes indicate the derivatives with respect to the variable £. Next we look for the solutions 
of these coupled total differential equations in r(£) and under some simplifying assumtions. 

After defining y = r 2 6' for the right hand side, eq.(13) can easily be recast in the form, y' = 
((v — w)/D)y, which admits a solution y = y exp[((v — w)/D)£\, or 

e' = (k/r 2 ).eM((v-w)/D)£], (14) 

where yo ( or k) is the integration constant. For simplicity, we concentrate here on the case when 
w = v , i.e., when the parameter w in the ansatz (11) becomes the convective velocity v of the system. 
This leads to & = k/r 2 . 
For the above choices, eq. (12) can be expressed as 

r " = K + \ Br + l Ar2 ' 
r 6 4 8 

which can be easily integrated to give 

(r') 2 = -% + -Br 2 + -Ar^ + h (15) 

r 2 4 4 4 



where (&i/4) is the constant of integration and B = 4(5 — a — Uq)/D, A = (8b/3D). Alternatively, 
by defining S = r 2 we write eq.(15) in terms of the variable S as 

(S') 2 = AVs E +BS 2 + k 1 S-4k 2 . (16) 

For the solitary wave solutions, we set k\ — k — 0, therby reducing eq.(16) to the form S' = 
S(A^S + B) l l 2 or eq. (15) to 



r' = (l/2)rVAr + B, (17) 
which can be integrated to give [18] r(£) as 



r(O = -(B/A).sec 2 (^-Z + Z ), 
for B < 0, and 



r(0 = (B/A).cosech 2 (^ + ^), 

for B > 0. Correspondingly, the solution of 0' = equation is taken as = constant (say 0q). 
Finally, the solutions of (4) in view of (11) becomes 

C(x, t) = -{B/A). exp(^ + St). sec 2 (^?e + Co), {B < 0), (18) 

and 



\/ B 

C(x, t) = {B/A). exp(^ + 5t).cosec/i 2 (^^ + f ), (-B > 0), (19) 

where ^o(= x o — vto) is the constant of integration and the same can be fixed from the initial 
conditions. 

B. Solution of eq.(5) 

Following the same procedure as for eq.(4) in the above subsection with ansatz (11) , the imaginary 
part of the resultant expression, in the present case, will yield the same equation as eq.(13). The 
real part however now becomes 

(v - w)r' + r5 = Dr" - Dr9' 2 + (a + U )r - br 3 . 
An equation analogous to eq.(16) for the present case can be derived as 

(S') 2 = AS 3 + BS 2 + k ± S - 4k 2 . (20) 



Under the same simplifying assumptions as made in the above subsection for the case of solitary 
wave solutions, the solution of 0' = equation is again taken as = 9 . The equation analogous to 
eq.(17) now becomes 
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r' = {l/2)rVAr 2 + B, (21) 
where B is the same as before but A = (2b/ D). Eq.(21) can be solved[18] to give 

r(0 = y/-B/A. sec(^£ + (B < 0), 

r(0 = ±y/B/A.cosech(^ + (5 > 0). 
Finally, the solitary wave solution of eq.(5) can be written as 

C(x, = sj-BjA. exp(i0 o + sec(^?e + Co), (5 < 0), (22) 

C(x, t) = ±\[b~/A. exp(i0 o + 5t).cosech(^-£ + Co), (5 > 0). (23) 
where £o is the constant of integration to be determined from the initial conditions. 

III. EIGENVALUE PROBLEM ASSOCIATED WITH EQ.(7) 
A. General results 

Since the Hamiltonian operator (8) is non-Hermitian, the eigenvalue A in (7) need not be real. 
Further for a symmetric potential, Ti is also not a PT-symmetric (a relaxed case of non-Hermiticity) 
one and hence may not admit [9] real eigenvalues. Naturally, the conjecture of Bender et al[9] for the 
reality of eigenvalues of PT-symmetric potentials (developed mainly in the context of the Schrodinger 
equation) is bound to show some limitations in this case. While the concept of pseudo-Hermiticity[ll] 
may be worth attempting for the present case, we have however been pursuing[12,13] an altogether 
different method to handle the complex Hamiltonian systems. From this point of view our approach 
is quite general and the concept of an extended complex phase space defined by (10) is used. In 
what follows we investigate the solution of an analogous D-R equation (7), in the sense that x and 
p in it are now complex. 

Note that V(x) in (7) in general could be complex just as x, p and ip are. Thus, we write 
V(x) = V r (x 1 ,p 2 ) + iV i (x 1 ,p 2 ),ip(x) = tp r (x 1 ,p 2 ) + iip i (x 1 ,p 2 ),X = \ r + iXi, and 

A. JL_iJL.Jl— 92 2 j 92 — 

dx dx dp 2 dx 2 dx\ dx\dp 2 dp 2 .' 
and use them in (7) to give 

B'[lpr,xixi '^'i ,rl Pr,xip2 VV,p2P2 Wpi,x\x\ ^^i,xip2 ^^Pi,P2P2\ ^[VV,xi ^VV,p2 ^Pi,x\ ^i,p^\ 

+ V r 1p r — Vilpi + iV r 1pi + iVilpr = \ r 1p r — Aj^j + i\lpr + i K^Pi -(24) 

Now we equate the real and imaginary part of this expression separately to zero. This leads to the 
following pair of coupled partial differential equation in ip r and ipi : 
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-D[-2^ Xlp2 + VVl*l - ^W 2 ] + V [-^r,p 2 + A,xA + V r 1pi + Vilpr = Xilpr + A r ^. 

Further use of the Cauchy-Riemann conditions for the analyticity of ip(x), viz., 

lftr,xi 1pi,p2 j VV,p2 ^i,x\i (25) 

leads to simpler forms of these equations, namely 

- 4.0^1*1 + 2vip riXl + V r ip r - Vitpi = X r tp r - Aj^i, (26) 

- AD%l) itXlXl + 2vip i>Xl + V r lpi + Vi1^ r = Xi^pr + X r tpi. (27) 
As for the case of Schrodinger equation[12], we make an ansatz here for the solution,viz., 

ip(x) = i) r + i^i = e 9(x) ;g(x) = g r (x 1 ,p 2 ) + ig i (x 1 ,p 2 ), (28) 

which gives ^ r (x 1 ,p 2 ) = e 9r cos gi,i/ji(xi,p 2 ) = e 9r sing, or g r = (1/2) ln(^ 2 + ip 2 ), gi = tan' 1 (ipi / 'ip r ) . 
Now, after using (28) in (26) and rationalizing the resultant expression with respect to the orthogonal 
functions cos(gj) and sin(gj), we obtain the following pair of coupled partial differential equations: 

9r,x lXl + {9r, Xl ) 2 - (gi, Xl ) 2 - (v/2D)g r , Xl + (l/AD)(X r - V r ) = 0, (29) 

gi,x ixi + 2g r , xl . gi , Xl - (v/2D)g hXl + (l/4£))(A i - V,) = 0. (30) 

Interestingly, the same set of equations is also arrived at if one rationalizes the eq.(27) using (28). 
Thus, for a given complex potential, eqs.(29) and (30) in which the original ansatz for ijj(x) is now 
transcribed into that for g(x) can be solved to obtain the real and imaginary parts of the eigenvalue 
A. Further, we shall consider the two situations- one when the parameter (s) of the potential V(x) are 
real or the other when they are complex. In the following we demonstrate the applications of these 
general results to the case of a simple harmonic oscillator corresponding to these two situations. 



B. Applications to complex harmonic oscillator potential 

We first consider the case of a complex harmonic oscillator potential, 

V(x) = ax 2 (31) 

where a is real. Using (10), we write V r (xi,p 2 ) = a(x 2 — p 2 ), Vi(xi,p 2 ) = 2ax\p 2 and make the 
following ansatz for g r and which is consistent with the Cauchy-Riemann conditions: 

g r (xi,p 2 ) = a 1 x 1 +a 2 p 2 + a 20 (x 2 1 -pl)+a n xip 2 ;g i (x 1 ,p 2 ) = -a 2 x 1 + a 1 p 2 -^a n ( y x 2 1 -pl) + 2a 20 x 1 p 2 , 

(32) 

where ctj, a^'s are real constants to be determined later. Substitution of (32) in (29) and subsequent 
rationalization of the resultant expression with respect to the powers of xi, p 2 and (xip 2 ) yields the 
following relations: 
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2a 20 + a\-a\- {v/2D)a 1 + (l/4D)A r = 0, (33) 

Aa\a 20 - 2a 2 a u - (v/D)a 20 = 0, (34) 

4a 2 o;2o + 2a±an — {v/2D)a u = 0, (35) 

Aa 20 an + Aa 20 a n = 0, (36) 

4« 2 2 - «n - (a/ AD) = 0, (37) 

«ii - 4 «2o + ( a / 4 -°) = 0. (38) 

Note that while eqs. (37) and (38) turn out to be identical, eq.(36) gives ana 20 = which implies 

that either a 20 = or an = (of course both of them can not be set equal to zero). Again a 20 = 



is ruled out as it leads to imaginary an from (37). Therefore, an = and a 20 = ±ya/lQD from 
(37). Eqs. (34), (35) and (33) immediately give a\ = (v/AD), a 2 = , and the real part of A as 

A r = T^VaD + (v 2 /AD). (39) 

Similarly, if one rationalizes (30) using (32) the same set of equations and hence the same values (as 
above) of the unknown parameters in the ansatz (32) are obtained except for the equation involving 
Aj. The imaginary part,Aj, of A now turns out to be zero, which implies only the real eigenvalue for 
the potential (31). Thus, for the potential (31), we have 

A r = T 2v^D + (v 2 /AD);\i = 0, (40) 



g r (x 1 ,p 2 ) = (v/AD)x 1 ± ^a/\QD(x\ - p\); gi(x 1: p 2 ) = (v/AD)p 2 ± ^a/ADx x p 2 
or equivalently, 



ijj(x) = exp[(v/AD)x ± ^a/lQDx 2 }. (41) 

Next, we consider the case when the parameter a of the potential (31) is complex, a = a r + iai. 
In that case we have 

V r {x u p 2 ) = a r (x\ -p\) - 2a i x 1 p 2 ; Vi(x u p 2 ) = a^xj - p\) - 2a r .x 1 p 2 . (42) 

As before, the use of ansatz (32) in (29) now yields the following set of equations after the rational- 
ization of the resultant exptression: 

2a 20 + a\ - a\ - (v/2D)a 1 + (l/4D)A r = 0, (43) 
4a 2 a 20 + 2a x a xl - (v/2D)a u = 0, (44) 
Aaia 20 - 2a 2 a u ~ (v/D)a 20 = 0, (45) 
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Aa 20 a u + 4a 2 ottn + (04 /2D) = 0, (46) 
4«2o - < ~ K/4£>) = 0, (47) 

«n - 4 «2o + K/4£>) = 0. (48) 

Eqs.(44)-(47) can be solved for the real parameters a±, a 2 , «2o and «u in the ansatz (32) and in 
terms of the complex parameter a of the potential (31) to give 

«i = (v/AD); a 2 = 0; a 20 = a+/4v / 2~D; a u = -a_/2v / 2~D, 
where a± = ^\a\ ± a r . Using these results in (43), one obtains A r as 

A r = T V2Da + + (v 2 /AD). (49) 

Similarly, the rationalization of eq.(30) using (32) and (42) will give the imaginary part, Aj, of the 
eigenvalue in a straightforward manner,viz., 

Xi = ±V2Da_. (50) 

Finally, the eigenfunction ip(x) in terms of g r (xi,p 2 ) and gi{x\,p 2 ) turns out to be 

g r {x 1 ,p 2 ) = (v/AD)x 1 + (a + /AV2D)(xl-p 2 2 ) - (a_/2v / 2D)x 1 p 2 ; 

9i{xi,p 2 ) = {v/AD)p 2 + (a-/AV2D)(xl-pl) + (a + /2V2D)x lP2 , 

or equivalently, 

tj}{x) = exp[(v/AD)x + (l/4v / 2D)(a+ + ia_)x 2 ]. (51) 

Before we highlight two special cases of these general results it can be seen that the eigenvalues 
for the non-Hermitian operator (8) are real as long as the potential parameter a is real. Once a 
becomes complex, then the complexicity of the eigenvalue may arise as a result of a i 7^ 0. 

Note that for the case when D = (1/2) and v = 0, above results for the complex harmonic 
oscillator trivially reduce to those derived by solving the Schrodinger equation in extended complex 
phase space (cf. Ref.(13), Sect. 3). Further, if we set = 0, a r = a = \a\ (or a_ = , a + = V2a), 
then it is not difficult to see that results (49)-(51) for the complex coupling reduce to that for the 
case of real coupling (cf. eqs.(40) and (41)). 

Mainly for the sake of a comparision we present here the results for the case of real harmonic 
oscillator, V(x) = ax 2 , obtained by solving (7) in a real phase plane. Again using the ansatz 
ip(x) = exp[g(x)], with g(x) = ax 2 + fix for the solution, it is not difficult to arrive at the following 
results for A and ip(x) : 

A = tVoD + (v 2 /AD), (52) 



ip(x) = exp[(v/2D)x ± ^Ja/ADx 2 }. (53) 

It can be seen that the results (52) and (53) differ slightly from (40) and (41) in terms of numerical 
factors. This is maily because of the requirement of analyticity of ip(x) used in the case of complex 
phase space. 
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IV. OTHER FEATURES OF THE DIFFUSION HAMILTONIAN (8) 

In this section we discuss some previously unnoticed mathematical features of the D-R Hamiltonian 
(8) (or its 'classical' version (9)). In the first part we derive an Ermakov system of equations for 
(8) and in the second, following the prescription advanced in our earlier work[15], we derive an 
equivalent integrable system of two, two-dimensional real Hamiltonians corresponding to the complex 
Hamiltonian (8). 

A. Ermakov system of equations and space invariant 

About 125 years ago, Ermakov for the first time demonstrated[17] the linkage between the solutions 
of certain type of differential equations via an integral invariant. This latter construct, now termed 
as 'Ermakov' (or 'Lewis') invariant in the context of classical mechanics, has played[19] an important 
role in the study of time-dependent systems and in the quantum[16] and other contexts[21,22] several 
new interpretations of this mathematical construct have been sought. Here, since the variable x 
characterizes the space dimension, we shall call this construct as the 'space invariant'. 

For the present purpose, we rewrite eq.(7) as 

V>"-7^' + g 2 (x)0 = O, (54) 

where 7 = (v/D), q 2 (x) = (X — V(x))/D , and look for its solutions in the form (ansatz) 

ip(x) = 4>(x)exp[if(x)]. (55) 

Now, after using (55) in (54) and equating the real and imaginary parts separately to zero in the 
resultant expression, we obtain 

0"-(f) 2 0-70' + g 2 = O, (56) 

/"0 + 2/V - 70/' = 0. (57) 

As before, eq.(57) after defining y = (f) 2 f, can easily be recast in the form y' = 73/, whose solution 
now becomes y = K.exp^x), or 

/' = ^/0 2 , (58) 
where k is the constant of integration. While the integration of (58) in the form 

f{x) = k f{e^'/<i> 2 {x'))dx\ (59) 

suggests a phase-amplitude connection in the present case, its use in (56) leads to a nonlinear 
differential equation, 

0" - 70' + g 2 = /tV 7 70 3 . (60) 

In order to derive the space invariant, multiply eq.(54) by and eq.(60) by -0 and subtract the 
latter to give, 

(00>" - 0» + 7(^0' - V>'0) = -/t 2 e 2rr 0/0 3 . (61) 

This expression, after using 2(0-0' — 0'-0) as the integrating factor, can immediately be integrated to 
give the space invariant, K, in the form 

K = K 2 (V/0) 2 + e - 2 ^(0V'-0V>) 2 . (62) 
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Note that the structure (62) satisfies (dK/dx) = 0, and hence termed as 'space invariant'. Clearly, 
it is a manifestation of the phase-amplitude connection (59). Although, K as such appears to be 
independent of the form of V(x) but the fact is that the role V(x) enters through (60) or for that 
matter via (59). The system of eqs. (54), (60) and (62) constitute an Ermakov system. 

B. An equivalent integrable system in two real dimensions 

In our earlier[15] work, we have exploited the analyticity property of the one-dimensional complex 
Hamiltonian H(x,p) by writing it in the form H(x,p) = H 1 (x 1 ,x 2 ,pi,p 2 ) + iH 2 (x 1 ,x 2 ,p 1 ,p 2 ) (cf. 
eq.(10)). This has resulted into a new class of two-dimensional integrable real Hamiltonian systems 
described by Hi(xi,x 2 ,pi,p 2 ) or by H 2 (xi,x 2 ,pi,p 2 ) for certain choices of the complex potential 
V(x). Here we shall demonstrate that this prescription also works for the complex Hamiltonian 
(9) and more so independently of the form of V(x) which is now expressed as V(x) = V r (xi,p 2 ) + 
iV i (x 1 ,p 2 ). 

Note that the 'classical' version(9) of the Hamiltonian (8) is complex even in the real phase plane. 
In the complex phase space (10), this version of Ti , expressed as H(x,p) = Hi(xi,x 2 ,pi,p 2 ) + 
iH 2 (xi,x 2 ,pi,p 2 ), becomes a function of two complex variables x and p and its real and imaginary 
parts turn out to be 

H 1 (x 1 ,x 2 ,p 1 ,p 2 ) = D(p\-x 2 2 ) -vx 2 + V r (x 1 ,p 2 ), (63) 

H 2 (x 1 ,x 2 ,p 1 ,p 2 ) = 2Dp 1 x 2 + vp 1 + V i (x 1 ,p 2 ). (64) 
Next, we compute the Poisson bracket [Hi,H 2 ]pb from 



r „ „, dH, dH 2 dH x dH 2 dH, dH 2 dH, dH 2 

H2]pB = ^7-^7 " + -dx~ 2 -W, ' d^"d^- (65) 

It can be noticed that for H 1 and H 2 given by (63) and (64) the poisson bracket (65) vanishes, 
if and only if the real and imaginary parts of V(x) satisfy the Cauchy-Riemann conditions, viz., 
V r ,x! = Vi,p 2 , ^»,xi — —Vr, P2 - In other words, this implies the analyticity of V(x). Further, vanishing 
of the Poisson bracket (65) also suggests that H\ and H 2 are in involution and independent in the 
sense[15] that vectors Vj = JV y Hj(y) for j = 1,2 and y = Xi,x 2 ,p±,p 2 turn out to be linearly 
independent for (63) and (64). Here J is the symplectic unit matrix. It may be mentioned that Hi 
and H 2 as given by (63) and (64) fulfil all the other requirements listed in Ref.(15) for a biharmonic 
function or for an auto-Backlund transformation and thereby enabling the integrability of the system. 

For the integrability of a two-dimensional, time independent system one expects the existence of 
another invariant in addition to the given Hamiltonian of the system. Such a second invariant, if 
becomes available, helps in quadrature. While the search for such an invariant has been there in 
the literature for many years and for many systems[19], the same has been constructed only for a 
few countable ones. In spite of many developments in terms of various methods for this purpose[19], 
there still remains scarcity of integrable systems. As a matter of fact, in some cases, even if the 
existence of the invariant is confirmed, its construction becomes a problem. From this point of view, 
the prescription suggested in Ref.(15) offers some kind of relief in the sense that at least for a certain 
class of Hamiltonian systems the second invariant becomes automatically available. In the present 
case, the vanishing of the Poisson bracket (65) not only suggest that H\ and H 2 are in involution but 
also that if Hi(xi, x 2 ,pi,p 2 ) is the Hamiltonian of the system, then H 2 (xi,x 2 ,pi,p 2 ) is the second 
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invariant of this system or the vice- versa, and thus implying the integrabilty of the system. No doubt 
the structures of the Hamiltonians, in this case, depart from the conventional ones but the fact is 
that such structures have also become desirable to understand some newly discovered phenomena in 
nature (see the citations in Ref.(15)). 

In view of the above general discussion, a few remarks about the Hamiltonian corresponding to 
the case of complex oscillator (cf. Sect. 3(B) ) are in order. Note that the potential V(x) in this case 
is an analytic function only when the parameter a is real. Corresponding to this case the equivalent 
integrable system in two real dimensions (cf. eqs. (63) and (64)) turns out to be 

H 1 (x 1 ,x 2 ,Pi,P2) = D(pl -x\) -vx 2 + a{x\ - pi), (66) 

H 2 (x 1 ,x 2 ,pi,P2) = 2DpiX 2 + vpi + 2ax x p 2 . (67) 

These forms of Hi and H 2 are worth comparable to any one of the examples studied in Ref.(15). 
As a matter of fact for the case of complex parameter the potenial V(x) is no longer an analytic 
function in the spirit of the above Cauchy-Riemann conditions, in stead it becomes the function of 
two complex variables a and x. For the analyticity of such a function one has to use a generalized 
version of Cauchy-Riemann conditions [23] which, in fact, involve the higher derivatives of V r and V^; 
whereas the computation of the Poisson bracket (65) requires only the first derivatives of V r and V{. 

V. CONCLUDING DISCUSSION 

From the point of view of learning more about the D-R equation, two slightly disconnected aspects 
of this equation are explored in this work. In the first part, we have studied the solitary wave solutions 
of the D-R equation with some specific types of quadratic and cubic nonlinearities, which, to the best 
of our knowledge, have not been investgated earlier in this context. With regard to the application of 
these results, it concerns the modelling part of the study of a nonlinear phenomenon. In fact, there 
appear now many situations in the fields of population biology or in different branches of physics and 
chemistry, where the results obtained here can be useful in offering the alternatives while accounting 
for the nonlinearity in such studies. 

We have restricted ourselves only to the study of solitary wave solutions of eqs. (4) and (5) by way 
of ignoring certain terms in eqs. (16) and (20). This is done mainly for simplicity ; otherwise one 
can as well retain all the terms in these equations and integrate the resultant expressions to obtain 
the traveling wave solutions of more general type. In fact, it is not difficult to arrive at the cnoidal 
wave solution in case of eq. (5) for certain choices of the parameters and in analogy with the KdV 
equation [20]. 

In spite of the fact that the Hamiltonian (8) is non- Hermit ian, it does not admit complex eigenval- 
ues in general (cf. Sect. 3 ). The complex eigenvalues can, however, be expected if the parameter(s) 
of the potential also becomes complex in addition to the phase space. In some sense this will bring 
the D-R equation closer to the Schrodinger equation in spite of the presence of the velocity-term in 
the former, particularly with reference to the large-t behaviour of their solutions. 

Two previously unexplored (to the best of our knowledge) aspects of the D-R Hamiltonian (8) 
are highlighted in the last section. These studies only hint to the richness of the mathematical 
content present already in the Hamiltonian (8) or in its 'classical' version (9). In analogy with other 
studies[16,19,21,22], a couple of possible interpretations of the constructed space invariant K (cf. 
eq.(62)) can be re-emphasized here in the present context. Firstly, since K is a space invariant and 
involves the solutions of both the eqs. (54) and (60), it can act as a geometric constraint with regard 
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to the validity of these solutions. Secondly, as the invariant K is the manifestation of a particular 
type of phase-amplitude connection (59) via the nonlinear eq. (60), its existence itself suggests[19,22] 
a nonlinear superposition principle in which the solution of (60) is expressible in terms of two linearly 
independent solutions of (54). Further applications of some of these results to the stellar structure 
studies are in progress. 
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